STM imaging of electronic waves on the surface of Bi2Te3: topologically protected 

surface states and hexagonal warping effects 
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Scanning tunneling spectroscopy studies on high-quality Bi2Te3 crystals exhibit perfect corre- 
spondence to ARPES data, hence enabling identification of different regimes measured in the local 
density of states (LDOS). Oscillations of LDOS near a step are analyzed. Within the main part of 
the surface band oscillations are strongly damped, supporting the hypothesis of topological protec- 
tion. At higher energies, as the surface band becomes concave, oscillations appear which disperse 
with a particular wave-vector that may result from an unconventional hexagonal warping term. 
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A new type of three-dimensional (3D) bulk insulating 
materials with surface Quantum Spin Hall Effect states 
protected by time reversal symmetry has been recently 
predicted [l|, and soon afterwards observed experimen- 
tally in BiSb bulk crystals [2]. Subsequently, Bi2Te3 
has been argued to be a similar three-dimensional topo- 
logical insulator (TI), exhibiting a bulk gap and a sin- 
gle, non-degenerate Dirac fermion band on the surface 
[3|. Indeed, recent angle resolved photoemission spec- 
troscopy (ARPES) confirmed that prediction |4]. In par- 
ticular, with appropriate hole-doping, the Fermi level 
could be tuned to intersect only the surface states, in- 
dicating fully gapped bulk states as is expected from a 
three-dimensional TI. While ARPES could confirm the 
nature of the band, it is still a challenge to demonstrate 
unambiguously the topologically "protected" nature of 
the surface state in Bi2Te3, or any other 3D TI system. 

In this paper we present scanning tunneling microscopy 
(STM) and spectroscopy (STS) studies on high-quahty 
doped Bi2Te3 crystals. First we show that the STS spec- 
tra exhibit remarkable correspondence to ARPES data, 
hence enabling us to identify each region of the local 
density of states (LDOS) measured. Second, by ana- 
lyzing the electron-waves (Friedel-oscillations) observed 
near cleavage steps, we show that within the main part of 
the surface state band oscillations are strongly damped, 
a hallmark of the strong supression of backscattering, 
hence supporting the hypothesis of a protected band. Fi- 
nally, we show that in the region in which the surface 
band is warped, pronounced oscillations appear, with a 
distinct nesting wave-vector. Possible influence of the 
bulk conduction band on the oscillations is also proposed. 

For the present study we used Sn and Cd doped sin- 
gle crystals of Bi2Te3 (see Fig. [T^ for crystal structure 
and Fig. [T]3,c for ARPES data). Nominal doping levels 
between and 0.27% for Sn, and up to 1% for Cd were 



incorporated to compensate n-type doping from vacancy 
and anti-site defects that are common in the Bi2Te3 sys- 
tem. Actual doping was determined separately using 
chemical and Hall-effect methods and were shown by 
ARPES 3 to be in excellent agreement with the rela- 
tive position of the Dirac point with respect to the Fermi 
energy. For example, undoped crystals exhibit a Dirac 
point at ^ —335 meV, 0.27% Sn doping yielded a Dirac 
point at ~ —300 meV, while for a typical ~ 1% Cd-doped 
crystal used for this study we observed the Dirac point 
at ~ —275 meV, equivalent to Sn doping of 0.45%. 

Samples were cleaved in vacuum of better than 5 x 
j^Q-io ^Q]-p^ g^jjfj quickly lowered to the ~9 K section 
of the microscope, where cryo-pumping ensures that the 
surface remains free from adsorbates for weeks or longer. 
Topography scans were taken at several bias voltages and 
setpoint currents (usually 200mV and lOOpA). A typi- 
cal spectrum, taken at a random position on the sur- 
face of a Cd-doped cleaved crystal is shown in Fig. [TJl. 
Many topography and spectroscopy maps were collected 
for both, the Sn-doped and Cd-doped samples with scan 
size as large as - (360 x 360) A^. Both were found to be 
similar to previously pubhshed STM studies of Bi2Te3 
and Bi2Se3 [1, @, [3|, including the typical triangular and 
clover-shaped defects. Weak atomic resolution could be 
obtained only near defects, presumably due to the overall 
high conductivity at the surface. 

To understand the various regimes of the LDOS curves 
of Bi2Te3, we compare STS spectra to ARPES data. 
Fig. Wp shows six constant-energy contours of the sur- 
face state band (SSB) of an undoped crystal. Starting 
from the "Dirac point" near E k, —335 meV, the con- 
tour changes from a circle to a hexagon, and becoming 
warped above about E « —100 meV. The inset to Fig.[T]: 
shows the full dispersion of the SSB in the two principal 
directions of the Brillouin zone (BZ), together with the 
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FIG. 1: DOS of Sn-doped Bi2Te3. a) Crystal structure show- 
ing three quintuple units, b) ARPES-measured Constant- 
energy contours of the SSB. Strength of DOS grows from blue 
(no DOS) to dark red (strong contribution to DOS). Arrows 
point to the position of maximum in DOS (note six fold sym- 
metry), c) Integrated DOS from ARPES. Here: Ep is the 
Fermi level; Ea is the bottom of the BCB; Eb is the point of 
the "opening up" of the SSB; Ec is the top of the BVB; and 
Eb is the Dirac point, d) Typical STS spectrum of a 0.27% 
Sn-doped Bi2Te3 sample similar to (c). Note the shift of all 
energies by ~ 40 meV as a result of the doping. 



portions of the bulk valence band (BVB) and the bulk 
conduction band (BCB) in its vicinity. Integrating the 
full i?(k) data Q over the BZ, we obtain the density of 
occupied states at the surface of the sample as a func- 
tion of energy. Fig. [H; is the result of this integration in 
the relevant energy range that includes the Dirac band 
which resides in the ~250 meV gap of this system. The 
directly measured LDOS of a similar crystal, as obtained 



by STM is shown in Fig. [TJi, where we could extend the 
data to positive bias. In both figures zero energy marks 
the Fermi level (Ep), Ea corresponds to the bottom of 
the BCB as measured by ARPES and Eb corresponds to 
the point where the SSB becomes warped (see Fig. [T]d). 
The linearly-dispersing Dirac band extends from its tip 
at Ee> to ~ Eb, while Ec denotes the top of the BVB. As 
can be seen from the ARPES dispersion in Fig. [1]:, the 
Dirac point is not exposed in this system, hence leaving 
us with only a relatively small range of pure linear LDOS 
expected from a Dirac dispersion. The remarkable agree- 
ment between the DOS obtained from ARPES and from 
STM suggests that no unusual tunneling matrix elements 
exist in this system to obscure the surface band. 




FIG. 2: a) Topography of the step. The total height of the 
step ~30 A corresponding to one unit cell. The height of re- 
gion II is ~20 A. Also shown is IQO-fold higher resolution to- 
pograph of region III where oscillations have been studied, b) 
Example of (raw data) oscillations at constant energy (-f 200 
meV) across the width of the step, c) Averaged spectroscopy 
as a function of distance from the step. The position of the 
step corresponds to "0" of the x-axis. Note the abrupt ter- 
mination of the LDOS oscillations at Eb- (Visible features 
near Ep are due to convolution with the electronic spectrum 
of the tip arising from its sharp geometry). 

As has been previously shown [1, 0, [l^], the elec- 
tronic properties of a surface state can be studied through 
the analysis of interference patterns of electron-waves 
(Friedel-oscillations) near defects. The simplest case is a 
step along a crystallographic axis, which in principle re- 
duces the interference problem to a ID phenomenon Q]. 
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Fig. [2^ is a topography, showing a step defect that was 
obtained in the process of cleaving the crystal. The step 
runs along the [100] direction, thus, oscillations will be 
searched primarily along the [120] direction. The down- 
step is marked as region-I in the figure. The total mea- 
sured thickness of the step between regions I and III is 
~30.5 A, in excellent agreement with the thickness of 
one unit cell [ll| (i.e. three quintuples). Examination 
of the topography also reveals a partial cleavage area, 
~6A thick, intermediate between regions I and III and 
we mark it as region-II. We will concentrate on the very 
flat region-Ill, in which this combined defect is expected 
to scatter quasiparticles and yield interference patterns 
that reflect the properties of the different bands that will 
be probed as a function of bias voltage. We mark the 
beginning of the step with a; = (we note that since the 
edge of the step is not atomically smooth, this choice is 
arbitrary to within a few A). We will explore the elec- 
tronic structure along a straight section of the step, a 
strip which is 80 A wide and 150 A long. 

Typical variation of LDOS oscillations across the step 
are shown for V = +200 mV in Fig. [2)3. Averaging simi- 
lar patterns over the width of the step for each energy in 
the interval between V = -400 mV and V = +200 mV, 
we obtain Fig. [2j:. Examination of this figure yields the 
following initial observations: i) Friedel-like oscillations 
[l3 ] that originate from the step are observed for all en- 
ergies above ~ Eb, which is the energy above which the 
SSB warps; ii) the period of the oscillations in that region 
increases with increasing bias voltage; and iii) the "nor- 
mal" oscillation pattern becomes strongly attenuated be- 
low ^ Eb with at most two visible peaks. 

We now turn to the analysis of the data. Fig. [3] 
shows a series of 8 energies ranging from high bias that 
includes the BCB, to near the Dirac point (which in- 
clude BVB). Panels a— e are in the regime of clear dis- 
persive oscillations. Here we fit the curves to a simple 
p{x) oc A{E)sin(2kx + (j))/x expression, appropriate for 
ID oscillations [8| of a surface band (with an amplitude 
A, and a scattering potential phase-shift (p). Attempts 
to fit the data with decays that are either much stronger 
(e.g. x~'^) or much weaker (e.g. x~'^'^) gave markedly 
poorer fits. In addition, pure bulk oscillations would re- 
quire an unobserved fc-divergence at Ea- Thus, the data 
strongly suggest that all the observed oscillations origi- 
nate from the SSB. While the oscillations seem to fit well 
the theoretical curve away from the step, the fit near the 
step is poor, especially close to x — 0. This is expected 
due to the step's height variation and its roughness [1] 
(Fig. [2^). In addition, the fit near the step may be dis- 
rupted due to a nondispersing peak at ~25 meV observed 
from below ~100 meV, down to the Dirac point. 

The amplitudes of the fitted oscillations are plotted in 
Fig. H^, and the respective values of k and (f> are plotted 
in Fig. [3)d. Also shown in Fig. Hja are the ARPES ex- 
tracted dispersion of the SSB for an undoped crystal in 
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FIG. 3: Sample averaged LDOS as a function of distance 
from the step for 8 energies. Red lines are fits to the data 
(see text), a), b), c), d), and e) correspond to energies above 
Eb showing well pronounced oscillations, f), g), and h) cor- 
respond to energies within the SSB, for which "h" may also 
reflect substantial influence of the BVB. 



the F-M and F-K directions. In the inset of Fig. IHd we 
use a —15 meV contour of the SSB to mark these two 
wavevectors as well as a "new" wavevector, knest which 
emerges above the energy in which the Fermi surface 
changes character from convex to concave. The disper- 
sion of knest agrees remarkably well with the dispersion 
extracted from the STM data below 70 meV, once the 
Dirac points are shifted from -335 meV in the undoped 
crystal to ^ —275 meV appropriate for the Cd-doped 
crystal. The associated phase changes by an amount of 
vr in that regime. Above 70 meV a break in the linear 
dispersion is observed, accompanied by a change of the 
phase back towards its initial value. We believe that at 
this energy the BCB influences the scattering, though the 
oscillations are still 2D in nature as we do not observe a 
change in the power of the decay. 

The abruptness with which the oscillations cease below 
~ —100 meV (i.e. below ~ Eb) strongly support the pre- 
diction that Bi2Te3 is a TI [sj. For any "normal" surface 
state band we would expect ~ a;^ ^-decaying oscillations 
to persist throughout the band. By contrast, SSB in TI 
systems should be protected from backscattering due to 
their definite chirality, and thus should not exhibit os- 
cillations, which are a consequence of the interference of 
incoming and backscattered waves. This "protection" re- 
sults from the definite chirality of the surface state. Thus, 
for scatterers that respect time reversal symmetry, any 
backscattering process should also include spin-fiip in or- 
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FIG. 4; Results of fits of oscillations (see text for details) for a 
Cd-doped sample with Dirac point at ~ —275 meV. a) Vari- 
ation of oscillations amplitude, b) Oscillation wave- vectors 
(full circles) and phase shifts (open circles). Triangles follow 
kr-M and kr-K, and diamonds follow knest, all three are ex- 
tracted from ARPES data of an undoped sample (see inset 
for example). Arrows point to the respective Dirac points. 



der to stay within the same chirality class. Since no such 
events are possible for non-magnetic defects, such as the 
step, no backscattering is allowed. However as is evident 
from Figs.[3f, g and h, the LDOS profile in the pure Dirac 
regime exhibits a nondispersing peak at ~ 25A and weak 
structure beyond it, suggesting at most, fast-decaying os- 
cillations below Eb and possibly a bound state near the 
step tl3|] . The idea of a bound state is also supported by 
the fact that the phase of the scattering potential changes 
by an order of ~ tt within the SSB. 

While the above analysis of the dispersion strongly sug- 
gests that the oscillations originate exclusively from the 
SSB, the BCB can still play a role in promoting oscilla- 
tions if we assume that due to its presence, quasiparticles 
can be scattered from k to -k via intermediate states in 
BCB. While calculations of such an effect are cumber- 
some and thus beyond the scope of this present work, it 
is not impossible that such processes maintain the 
type oscillations Since processes involving BCB are 
expected to affect the dispersion of the oscillations, it is 
reasonable to conclude that above ~ 70 meV the observed 
change in dispersion is due to BCB effects ( Fig. Hb)- 

Maybe the most surprising observation is the emer- 
gence of linearly dispersing oscillations between ^ —80 
meV and ~ 70 meV within the surface-state Dirac band. 
Here we propose that this is a direct consequence of the 
warping of the SSB, as its contour changes character from 
a fully convex shape (smoothly evolving from the Dirac 
point to a circle, and to a hexagon higher up in the band) 



to a concave hexagram (Fig. [T]d). In addition, above 
~ Eb DOS shift from the tips of the hexagons to the 
warped sides of the hexagram allowing for a new, un- 
protected, nesting wave-vector to appear. The inset in 
Fig. depicts this nesting wave-vector in the ARPES 
data. The linear dispersion extracted from the STM 
data, in perfect correspondence to knest extracted from 
ARPES, suggest that this is the only possible wave- vector 
for oscillations within the warped SSB. This is easily un- 
derstood if we assign a small z-direction component for 
the spin (Sz)- Such a component will have to alternate 
in sign between warped hexagonal sides, hence the only 
possible nesting wave-vector that respects Sz is the ob- 
served knest- Thus, the above analysis strongly suggests 
that the oscillations that emerge from the otherwise "pro- 
tected" SSB are a result of an unconventional hexagonal 
warping term [l5[ of the Bi2Te3 TI - system . 

Note added: Following posting of the first version of 
this paper on the archive, Liang Fu posted a paper [isj 
explaining the correlation between LDOS oscillations and 
convexity of constant energy contours. The revised ver- 
sion of this paper emphasizes this observation. 

Discussions with Qin Liu, Xiaoliang Qi, Shoucheng 
Zhang, and especially Liang Fu are greatly appreciated. 
This work was supported by the Department of Energy 
Grant DE-AC02-76SF00515. ZA is supported by Stan- 
ford Graduate Fellowship. 
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